Quantum order by disorder in frustrated diamond lattice antiferromagnets 
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We present a quantum theory of frustrated diamond lattice antiferromagnets. Considering quantum fluctua- 
tions as the predominant mechanism relieving spin frustration, we find a rich phase diagram comprising of six 
phases with coplanar spiral ordering in addition to the Neel phase. By computing the specific heat of these or- 
dered phases, we obtain a remarkable agreement between (/c, /c, 0)-spiral ordering and the experimental specific 
heat data for the diamond lattice spinel compounds MnSc2S4, C03O4 and CoRh204, i.e. specific heat data is 
a strong evidence for (/c, /c, 0) -spiral ordering in all of these materials. This prediction can be tested in future 
neutron scattering experiments on C03O4 and CoRh204, and is consistent with existing neutron scattering data 
on MnSc2S4. Based on this agreement we infer a monotonically increasing relationship between frustration and 
the strength of quantum fluctuations. 



Introduction. In insulating magnetic materials, new phases 
of matter may be found by letting local exchange interactions 
compete. In such situations, the spins are said to be frus- 
trated and intriguing new phases such as ordered phases with 
coplanar or spiral ordering or "spin liquid" paramagnets can 
arise. HI] In any given frustrated material, the ground state may 
be determined by identifying the primary mechanism reliev- 
ing the frustration. While extrinsic mechanisms, such as small 
dipole interactions, disorder or lattice distortions [2], may be 
important, perhaps the most interesting possibility is when 
temperature or quantum fluctuations alone relieve the frustra- 
tion, a process termed "order by disorder". [3] 

In this light, recent experiments which unveil strong frustra- 
tion in spinel compounds AB2X4, with magnetic ions occupy- 
ing the A- sites, are particularly interesting. Here the A- sites 
form a diamond lattice of spin S = f , 2, | local moments. 
Important examples include seven diamond spinels (|4] and 
references within), four that order (MnSc2S4, MnAl204, 
C03O4, CoRh204) and three that do not down to the lowest 
temperatures studied (C0AI2O4, FeAl204, FeSc2S4). When 
the moments order, the ordering temperature, Tc, is low com- 
pared to the Curie- Weiss temperature, 6cw, with frustration 
parameters [5], / = varying from 1.2 to 10. The mag- 

netic ordering in one of the ordered materials, MnSc2S4, has 
been identified as an exotic (/c, /c, 0) coplanar spiral via ex- 
tensive neutron scattering experiments [4] while the magnetic 
ordering patterns of other ordered diamond spinels are not de- 
termined yet. Given that a diamond lattice is bipartite, this 
ubiquitous evidence for frustration is highly unexpected. 

In combination with their frustrated magnetic properties, 
diamond spinels also have unusual temperature dependence 
of specific heat. Remarkably, among the four materials that 
order, their specific heat data share the same unusual behav- 
ior below Tc (see Ref.[4]). Instead of a pure power-law 
expected for incommensurate magnetic ordering in three di- 
mensions, two inflection points are observed. The three that 
do not order also share the same characteristic specific heat, 
but is quite different from those that order. These materials 
display a T^-^ power-law over a decade in temperature. [4] 

Following these experimental discoveries, the classical 
Heisenberg model on the diamond lattice with the nearest 



(nn) and next-nearest neighbor (nnn) exchange interactions 
has been studiedl^. It was demonstrated that the frustra- 
tion arises from the nnn interactions that couple spins within 
each of the two face-centered cubic (fee) sublattices of the 
diamond lattice structure (6*]. This coupling creates a highly 
degenerate set of classical coplanar spirals whose propagation 
vectors form a continuous surface in momentum space. Re- 
lieving this classical ground state degeneracy by thermal fluc- 
tuations was then found to produce a rich phase diagram at 
the classical level|l6|], including the (/c, /c, 0) spiral phase dis- 
covered in the neutron scattering experiments on MnSc2S4. 
This put diamond spinels in a promising class of materials in 
which (thermal) "order by disorder" may be experimentally 
observed. However, this classical picture may not be suf- 
ficient to describe possible effect of quantum fluctuations in 
these materials with relatively small spin 6* = |, 2, | and the 
specific heat data below Tc. 

In this letter, we present a quantum theory of frustrated dia- 
mond lattice antiferromagnets. We find that quantum fluctua- 
tions act as an order by disorder mechanism to produce a sim- 
ilar but richer phase diagram compared with that obtained ex- 
clusively from thermal fluctuations. In particular, focusing on 
the ordered states, we demonstrate that the characteristic sig- 
natures in the specific heat data of three of the ordered materi- 
als (MnSc2S4, C03O4 and CoRh204) can be explained if the 
magnetic ordering pattern in these materials is (/c, /c, 0) -spiral 
selected by strong quantum fluctuations. Thus, we argue that 
specific heat data is a convincing evidence for (/c, /c, 0) -spiral 
ordering in all these materials. Based on this comparison, we 
show how frustration and quantum fluctuations are strongly 
intertwined in these systems. Finally, we discuss the implica- 
tions of our results on future neutron scattering experiments 
on these materials. 

Keeping in mind that frustration arises from the nnn ex- 
change interactions!^, we begin with the quantum Heisen- 
berg model: 



H 



J2 



(1) 



where are spin-iS operators at site i, Ji > is the exchange 
coupling on the nn links (between sites on different fee sub- 



lattices) and J2 > is the exchange coupHng on the nnn Hnks 
(between sites on the same fee sublattice). By studying the 
large- A/^ limit of the Sp(A^) generalization of this modelj?!], 
we study the role of quantum fluctuations as a controlled ex- 
pansion in . The advantage of this method is that, unlike 
the large-iS expansion|8], the results are non-perturbative in 
the spin magnitude S (strength of quantum fluctuations) and 
the coupling constants (Ji and J2). The resulting phase di- 
agram, which indeed exhibits order by disorder, is presented 
in Figdl For ^ > 1/8, the energy of the degenerate set of 
classical states is given by the grey-to-black pixellated surface 
on the right hand side of the figure. Zero point energy correc- 
tions due to quantum fluctuations lift this degeneracy (grey- 
to-black pixels represent higher-to-lower energy spin config- 
urations) and only the black pixels with the smallest energy 
remain degenerate. The phase diagram was obtained based on 
the resulting magnetic order of the selected ground state. 

Fig [T] also shows that quantum and thermal fluctuations lift 
the classical ground state degeneracy differently. While ther- 
mal fluctuations entropically select an ordering pattern, quan- 
tum fluctuations select an ordering configuration purely from 
energetic considerations and so need not select the same state. 
Comparing our results with a study where only entropic ef- 
fects were considered fl, we notice that both kinds of fluctua- 
tions select states along similar high symmetry directions such 
as (/c, /c, k), (/c, /c, 0) and (/c, 0, 0). However, similar states are 
not always present in the same range of . In addition, within 
numerical accuracy, quantum fluctuations do not always lift 
the degeneracy, or lift it only partly as in the "degenerate", 
"circle" and "cross" phases. This is in contrast to thermal fluc- 
tuations, where points of lowest energy are always selected 
when entropic effects are considered. 

Focusing on the diamond spinels that magnetically order at 
low temperatures, we find remarkable agreement between the 
measured heat capacity and our large- A/" theory in a phase with 
(/c, /c, 0) ordering. As shown in Fig 121 ^ for CoRh204, C03O4 
and MnSc2S4, have two characteristic inflection points before 
reaching Tc, a feature that is best reproduced by the (/c, /c, 0) 
spiral ordering in the large- limit. Consequently, we pro- 
pose that all of these spinels most likely have (/c, /c, 0) spiral 
ordering. This result is in agreement with neutron scattering 
experiments on MnSc2S4.|4] Future neutron scattering exper- 
iments on CoRh204 and C03O4 could verify this prediction. 

Through the remarkable fitting of the large- A/^ theory to 
the heat capacity data, we also find a simple relationship be- 
tween the empirical frustration parameter, /, and the strength 
of quantum fluctuations given by the large- A^ effective spin 
length parameter = ^jf - (here is held fixed in the A/" ^ 00 
limit). To describe adequatly the experimental specific heat 
data of the more frustrated (larger /) compounds, we need 
to include stronger quantum fluctuations than required to de- 
scribe the moderately frustrated ones. As shown in Fig 121 the 
spin I system MnSc2S4 with frustration parameter / = 10 is 
best fitted using = 1.55 while the spin | system CoRh204 
with / = 1.2 is best fitted with Hi = 2.5. Thus for the diamond 
spinels that order at low temperatures, quantum fluctuations 



19 

CD 



3 



A ^ • # • I 



(k,0,0) 



4 points 
around 
(k,0,0) ^ 
00 



(k,k,0) 



CM 



6 points 
around ■- 
(k,k,k) 



(k,k,k) - 



Neel 



(k,0,0) 



Cross at / T -«JfM^ 
(k,0,0) / ^^'"1*; 

. £ points 
around 
(k,k,k) 
[include 
(k,k,0)] 

Circle 
around 
(k,k,k) 




(k,k,k) 



Degen^ate 



Neel 




ek 

k-space 



271 



FIG. 1 : Structure of the ground states as a function of ^ after the in- 
clusion of quantum fluctuations. The 0^ versus Ok plots on the right 
are planar views of the somewhat distorted and/or punctured sphere 
forming the set of degenerate classical ground states. 0^ and Ok are 
the polar and azimuthal angles spanning all directions in momentum 
space. Darker points represent lower energy. On each plot, red dots 
or lines give examples of states selected by quantum fluctuations. To 
contrast the effects of quantum and thermal fluctuations, states se- 
lected by entropic effects are summarized on the left. 1 6] 



correlate with frustration much more than with the physical 
spin representation. 

General Sp(N) mean-field free energy. The Hamiltonian 
describing the interactions between spins on the diamond lat- 
tice is given by Eq.©. The SU(2) spin symmetry of the 
Hamiltonian is generalized to Sp(A^) by first recasting the 
spins using the bosonic representation Si = ^bl^aafdbip 
where = {T, i} labels two possible spin states of each 
boson and then by introducing N flavors of such bosons on 
each site. In order to keep the physical Hilbert space of spins. 
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FIG. 2: Comparison of specific heat data of CoRh204, C03O4, 
MnSc2S4 and the theoretical large- specific heat of the (/c, /c, 0) 
spiral ordering with J2/J1 = 0.6. Here 1/k, = A^/25'eff, held 
fixed in the large- limit, gives the strength of quantum fluctu- 
ations and increases monotonically with the frustration parameter 
/ = IBcwI/Tc. The inset shows the subtraction of the nuclear 
contribution with a constant "background" from the specific heat of 
MnSc2S4 using Ci/T = A^T^e^^^^^ / (1 + e^^/^))^ + Co. 



a constraint on the number of bosons given by = bj^b'^ = 
25'eff = hzN where m = 1, must be imposed at each 
site. Note that N = 1 corresponds to the physical limit Sp(l) 
= SU(2). The action of the SpCA/") generalized model is then 
given by 



Jo 



where Aij 



^A,jA,j^X,{bZbTa-nM)} 
'^Ta^% (^cxf3 ^mm' IS the S^{N) gener- 



alized antisymmetric tensor of SU(2)) and the chemical poten- 
tial \i keeps the average number of bosons fixed to = nN 
at every site. The mean-field action is then obtained by decou- 
pling the quartic boson interaction in S using the Hubbard- 
Stratonovitch fields Qij = —Qji directed along the lattice 
links so that one obtains Qij = {Aij) /N at the saddle point. 
The mean field solution becomes exact in the large- A/^ limit 
where A/" ^ 00 is taken while n = rib /N is fixed. We also 

introduce the paramerization 6^ = ^/N Xia b^ ) where 
rh = 2, A^ to allow for the possiblity of long-range order 
that occurs when Xia 7^ 0. Consequently, after integrating 
over the bosons, and reseating Qij and A by hz, Xia by ^/K, and 
the temperature by hz'^ , we obtain the mean-field free energy 

^ E 1' - Qvi^c.0x:^^*p) + C.C.) 



where /eff 



(- + l))+/eff 



2kBT\n{l - e 



(2) 



"), and 



^(Q, A) are the eigenvalues of the mean-field Hamiltonian. 



Note that the chemical potential is now taken to be uniform 
since each site has the same number of nn and nnn links. In 
general, magnetic ordering x^q, 7^ occurs in the semiclassi- 
cal limit at larger k while quantum paramagnetic phases are 
obtained when k is small. 

Classical ground state. In the classical limit hz ^ 00 at 
T ^ 0, one can show that Q^j = Capx] 
classical energy is given by 



Xjp, so that the 



EJij_ 
2 



^af3Xia^jf3 



J'-l)- (3) 



Minimizing Eq.© with respect to and A^ is equiv- 
alent to determining the classical ground states of Eq.O 
provided the solution has = 1. Rewriting in 

terms of a quadratic form in the classical unit spin vectors 
— transforming to k- space unveils two 

bands e±(k) = J2(A^(k) - 1) ± ^ JiA(k) where A^{k) = 

4{nn=x,^,. cos2 ^ + Uu=x,y,z ^in' The minimum 

eigenvalue is obtained in the lower band e_ and is unique at 
k = for < I but highly degenerate, corresponding to a 

surface in k-space, for > |. For | < -j^ < |, the space 
of ground states resembles the surface of a slightly deformed 
sphere. For larger j^, the size of the "sphere" increases and 
eight holes centered around the (k^k^k) directions begin to 
puncture its surface (see also Ref.[6]). The solution of the 
classical limit is then completed by finding for each degen- 
erate eigenstate, labeled by its k value, its real space spin 
configuration and chemical potential. In terms of our spinor 



representation, we obtain 



and 



= _l^gT2(k•r^+2^(k)) ^hej^e A and B label the two 
fee sublattices and 6>(k) = arctan(tan ^ tan ^ tan ^). 



Varying E^ with respect to xl 

\c - 1 
^0 — z:^ 



gives the chemical potential 



Effect of quantum fluctuations. Given the above solution 
to the classical ground states, consider expanding the ground 



state energy E in powers of 1 / /t: so that E 
This leads to the quantum correction 



E^ 



= ^uj^{Q',X')-X'Ns 



where Ns is the number of lattice sites and Q^, the classi- 
cal values for the link variables, are given by Q^f^^^^ = 



sin( 



k-(rj-rO 



) and Q^,f/^^ 



± cos( 



k-(r/-rO±^(k) 



Since is only dependent on the difference between two 
sites, we can Fourier transform back to k-space and solve for 
E'^ analytically as a sum over wave- vectors. In practice, this 
energy correction is computed using an adaptive Monte Carlo 
integration method and the new ground state is found by sam- 
pling the surface of equal energy. The first order quantum 
corrections dramatically alter the topology of the degenerate 
ground state manifold reducing the k-space surface of lowest 
energy to only points or lines represented by the black pixels 
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inFigH] For 0.125 < 4^ < 0.18, the "sphere" of equal energy 
remains surprisingly degenerate; for 0.18 < ^ < 0.25, the 
eight {k^k^k) directions are selected; for 0.25 < ^ < 0.5, 
states labeled by k points forming eight circles around the 
(/c, /c, k) directions are chosen; for 0.5 < -j^ < 0.65, each 
circle gives way to six points around each {k^k^k) direc- 
tion (among which three points are (/c, /c, 0) directions); for 
0.65 < -j^ < 0.75, the states labeled by k points form a 
degenerate cross centered around each (/c, 0, 0) direction and 
finally for j^. > 0.75 states labeled by points along the six 
(/c, 0, 0) directions become the states of lowest energy. 

Specific heat and comparison to experiments. To compare 
the above order by disorder predictions with current experi- 
ments, we compute the specific heat. At finite temperatures, 
we assume that the phase and amplitude of both bonds and 
condensates can vary but that these changes are spatially uni- 
form. We write Qij = R e^^ and Xia, = \fr e*^ 
where and are bond and condensate values in one of 
the T = spin configurations chosen by quantum fluctua- 
tions. For a given T and effective spin length, R{Q^ A, T), 
r((5, A, T), ^ and C are obtained from self-consistent saddle 
point equations. For spiral configurations (k^k^k), {k^k^O) 
and (/c, 0, 0), the ground state is magnetically ordered at T = 
if > where, for example, tZc = 0.17, 0.094, 0.087 for 
^ = 0.2,0.6,0.85. 

Using Qij and Xia in Eq.©, we obtain the specific heat 
C = — T|^. In the limit of very low temperatures, R and r 
are approximately T-independent and C ~ as expected for 
3D antiferromagnets. However, as temperature is increased, 
R and r become T-dependent and the specific heat departs 
from its behavior. As shown in Fig|2l the specific heat ob- 
tained from the (/c, /c, 0) spiral ordering with = 0.6 agrees 
well with the experimental specific heat data of CoRh204 and 
C03O4 and the fit looks reasonable for MnSc2S4. For these 
three materials, ^ presents two inflection points before Tc, a 
feature that is best reproduced by (/c, /c, 0) ordering and is ro- 
bust for a finite range of ^ as long as the magnetic ordering 
remains the same. Other ordering wavevectors do not repro- 
duce as nicely the characteristic temperature dependence of 
the specific heat data. The two inflection points arise due to 
the temperature dependence of the magnon spectrum. The 
first point arises as the spin-wave velocity departs from its 
T = value while the second point appears when the tem- 
perature becomes of the order of the magnon bandwidth. Re- 
garding MnSc2S4, the low temperature part of the experimen- 
tal specific heat has substantial nuclear spin contributions that 



we subtract using a simple two-level system formula as de- 
scribed in the figure caption of Fig|2l Also, Ref.[6] pointed 
out that the neutron scattering experiment on this material sug- 
gests =0.85 and implies that a third neighbor coupling J3 
is necessary to stabilize (/c, /c,0) order for this value of 
We expect that the latter and the error arising from the sub- 
traction of the nuclear contribution are the reasons behind the 
less satisfactory fit for MnSc2S4. 

Conclusion. We presented a theory of frustrated diamond 
lattice quantum antiferromagnets. Considering quantum fluc- 
tuations as the predominant mechanism relieving spin frustra- 
tion in this spin system we found a rich phase diagram con- 
sisting of six phases with coplanar spiral ordering in addition 
to the Neel phase. By comparing specific heat curves found 
in the large- A/" mean-field theory with data obtained from 
CoRh204, C03O4 and MnSc2S4, we propose that they all 
share the same magnetic order in their ground state: a coplanar 
spiral with propagation vector (/c, /c, 0). Note that the neutron 
scattering data is currently available only for MnSc2S4. From 
the fit in Fig|2l we conclude that a remarkable correlation ex- 
ists between the strength of quantum fluctuations, measured 
by the effective spin magnitude k, = and the empirical 

frustration parameter / = while hz and the physical spin 
magnitude S = | , | of the magnetic ions appear only loosely 
related. We expect future neutron scattering experiments will 
verify our predictions. 
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